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Abstract Atom gravimeters are equivalent to non-multi- 
level corner-cube gravimeters in translating the gravity 
signal into the measurement result. This enables de- 
scription of atom gravimeters as LTI systems. The sys- 
tem's impulse responses by acceleration, velocity, and 
displacement are found to have the shape of triangle, 
meander, and the Dirac comb resp. The effects of in- 
homogeneous gravity field are studied for constant and 
linear vertical gradients and self-attraction of the in- 
strument. For the constant gradient the effective mea- 
surement height is below the top of the trajectory at 1/6 
and 7/24 of its length for the fountain and the release 
types of the instruments resp. The analysis is expanded 
to the gravimeters implementing the Bloch oscillations 
at the apex of the trajectory. In filtering the vibrations 
these instruments are equivalent to the first-order low- 
pass filters, while other atom gravimeters are equivalent 
to the second-order low-pass filters. 

Keywords Atom gravimeter • LTI system • Impulse 
response function • Vertical gravity gradient 

1 Introduction 

In 1991 Mark Kasevich noticed that the phase shift of 
the atom interferometer is quadratically proportional 
to the free-fall time of the atoms [TJ [2] , thus realiz- 
ing the first gravimeter with cold atoms as test mass. 
The following years saw rapid progress in this type 
of instruments, which now approach in accuracy and 
even exceed in sensitivity the best corner-cube gravime- 
ters [3n6l. The atom and the corner-cube instruments 
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measure gravity acceleration in significantly different 
ways, with almost non-overlapping sets of systematic ef- 
fects, so bringing both types of the instruments together 
for comparisons is very beneficial for metrology. Similar 
to corner-cube instruments, many systematic effects of 
atom gravimeters can be estimated only by modelling 
and computer simulation. To validate the corrections 
for the effects it's therefore important to compare re- 
sults obtained with different models. In this paper we 
outline atom gravimeters from the viewpoint of the the- 
ory of linear systems. In particular, we find the impulse 
response functions of atom gravimeters and apply them 
to the analysis of the effects of the inhomogeneous grav- 
ity field. 

2 Atom gravimeter as LTI system 

In atom gravimeters the positions of the free falling 
atoms are related to the phases of the light field gener- 
ated by two counter-propagating lasers. In the 3-pulse 
gravimeters the phases fa, fa, fa of the field correspond 
to the positions of the atoms at the moments t\, £2, £3 
separated by the time interval T. The interferometer 
output in these gravimeters is [7|[8j 

Acj> = kgT 2 + fa - 2fa + fa, (1) 

where k is the wave vectoiQof the light field. The quadratic 
term of ([I]) is caused by the Doppler effect, as the ac- 
celerating atoms see the static light wave as linearly 
increasing frequency. The term can be cancelled by con- 
trolling the frequencies of the lasers during the free fall. 

1 In later publications this parameter is often called the 
effective wave vector k e ({ to highlight the cumulative action 
of two counter-propagating beams. We use simplified notation 
to avoid allusions to the effective measurement height h e {{ 
considered further in this paper. 
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The measurement of the acceleration in atom gravime- 
ters consists in the experimental search of the frequency 
increase rate a that compensates the Doppler shift. The 
acceleration is obtained from the a as 19 



g = 2ira/k. 



(2) 



The right a makes A(f> equal to for any T [fjj , so from 
([I]) it also follows 



fJ 



2^2 + fo)/kT 2 



(3) 



As z = 4>/k, the same acceleration could also be found 
from the coordinates, if they were available [101: 



g = (z 3 -2z 2 + Zl )/T 2 



(4) 



Substituting the distances traveled by the atoms from 
the coordinate z\\ S\ = z 2 — Z\, S2 — z 3 — z\ and the 
corresponding time intervals T\ — T, T 2 — 2T into Q 
reveals the equivalence of the 3-pulse atom gravimeters 
with the 3-level measurement schemas (hg[3]) used in 
some corner-cube instruments, e.g. [Ill: 



9 = 



S2 
T 2 



Si 

T1JT2- T\ 



(5) 



Let the acceleration of atoms during the free fall change 
like g(t). The distance relates to the acceleration via 
double integration, so ^ can be put as 

- / Io T Io9(r)drdt £ fJ g(r)drdt \ 2 
2T T J T' 



9 = 



(6) 



The above formula relates the acceleration of atoms g(t) 
to the measured acceleration g linearly, i.e. linear com- 
bination of partial accelerations 



g(t) = a gi (t) + bg 2 (t), 
translates into 
9 = ag~i + bg 2 , 



(7) 



(8) 



where g 1 and g 2 are the results of independent measure- 
ments of gi (t) and g 2 (t) . This linearity plus the time in- 
variance of ^ (meaning that T\ and T 2 are pre-dehned 
for a drop) enables treatment of atom gravimeter as lin- 
ear time-invariant (LTI) system . In such a system 12 
the input and the output signals are connected by the 
convolution operation: 



+00 



g{r)h{t - r)dr, 



(9) 



where h(t) is the impulse response function of the sys- 
tem. Applied to absolute gravimeters, the convolution 
enables presenting the measurement result in the form 



r 2T 



9 = 



g(t)w g (t)dt, 



(10) 
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Fig. 1 The acceleration impulse response of atom 
gravimeter as reaction on the Dirac 5-function: 
a - S(t) as input signal; 

b - the Heaviside unit step function u(t) — j 5(t)dt; 
c - the unit ramp function r(t) = / / 5(t)dtdt; 
d - the weighting function of a double integrator as 
r(t) turned backward w.r.t. the end of integration; 
e - the weighting function of atom gravimeter com- 
bined of two scaled and shifted double integrators. 
The process is illustrated in more details in [13]. 



i.e. as the acceleration of the atoms weighted over the 
measurement interval. Here we have g — g(2T) and 
w g (t) = h(2T—t), as the measured gravity is attributed 
to the end of the measurement interval. Therefore, the 
gravimeter's weighting function is its impulse response 
function turned backwards with respect to the end of 
the measurement interval. The impulse response can 
be found as reaction of the system on the Dirac delta 
function S(t) [12 . As ^ has two double integrators, 
we first find the reaction of a double integrator on S(t). 
The first integral of S(t) is the Heaviside step function 
u(t), the second integral is the unit ramp function r(x) 



"(*) = 



t < 0, 
t t > 0. 



(11) 



The equivalence of the double integration to the con- 



volution with the unit ramp (111 follows from the inte- 
gration by parts rule f Q udv = [uv]^ — J Q vdu : 



j [ g(r)dr dt 

Jo Jo v ~y' 



(12) 







t f g(r) dr 




. Jo 


Jo 



tg(t)dt= j g(t){T-t)di. 
Jo 



Two ramp functions combined according to the formula 
Q yield the the gravimeter's weighting function by ac- 
celeration w g (t): 



Wg(t) 



r(2T-t) 
2T 



I 

T 2 

2^ t_ 

T T 2 



r(T-t) 



< t < T, 
T < t < 2T. 



(13) 



Some systematic effects are more convenient to analyze 
in terms of the atoms' velocity or displacement rather 
than acceleration. Applying again integration by parts 
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explaining why the constant acceleration translates un- 
changed into the result, while the constant velocity or 
displacement yield no measured gravity. 



3 Analysis of systematic effects using response 
functions 

If the atoms' acceleration changes during the free fall 
like 



At) 



.(*) 



h 



im. 

T 2 



S(t-2T) 
T 2 



1 



I i 



<S(t) S(t-T-2T') 



g(t)=g + Ag(t), 



(19) 



T(T+T') T(T+T') 

T T' Tl t 



2T 



1 I 



the additional component caused by the term Ag(t) can 
be found as 



Sjt-T) 



S(t-T) S(t-T-T') 
' T(T + T') T(T+T') 



/■2T 

A^= / Ag(t)w g (t)dt, 
Jo 



(20) 



a b 
Fig. 2 Weighting functions of atom gravimeters (top to bot- 
tom): by acceleration, by velocity, by displacement, 
a - 3-pulse atom gravimeter 

b - atom gravimeter implementing Bloch oscillations 



to (10), we get 

g= / g(t)w g (t)dt 



where w g (t) is the weighting function of the gravime- 
ter by acceleration. If the disturbance is expressed in 
terms of velocity or displacement, the additional com- 
ponent can be found using the corresponding weighting 
functions w v (t) or w z (t) as 

Ag^ = T- 2 (^J AV{t)dt- AV(t)dtj , (21) 



(14) 



w g {t) / g(t)dt 



2T 



2T 



dWgjt) 

dt 



g(t)dtdt. 



Ag = T 



r 2T 

2 I Az 



z(t) (S(t) - 2S(t - T) + S(t - 2T)) 



As w g (0) = w g (2T) = 0, and the velocity V(t) 
J g(t)dt, the above expression simplifies to 



= ( Az(0) - 2Az{T) + Az(2T) ) T 



dt 
(22) 



2T 

n = / V{t)w v {t)dtdt, 
h 



where 



w v (t) 



dw g{t) 
dt 



(15) 



(16) 



It's interesting to observe that the formulas (22 1 and 
Q are alike, as the gravimeter's weighting function by 
displacement is just the Dirac comb sampling the con- 
tinuously changing coordinate in three points. 
If the disturbance is expanded like 



N 



is the weighting function by velocitj]^] The weighting 
function by displacement can be found likewise, leading 
to 



Ag(t) = '£ i a n t n , 



(23) 



its influence on the measured gravity, according to ( 10 1, 
can be found as 



, \ d . . d 2 



(17) 



These functions are shown on the fig(2^i. The following 
properties hold true 

2T p2T p2T 

w g {t)dt=l, w v (t)dt = 0, w z (t)dt = 0, (18) 











2 The function similar to w v (t) is known in atom interfer- 
ometry as sensitivity function 14 . 



N 

Ag = y^a„C», 

i—n 

where C n is the n-th moment of the w g (t): 

r2T 2 »+2 - 2 

C " = /o ^^= ( n + l)(n + 2) r " 



(24) 



(25) 



Similar formulas can be derived also for w v (t) and w z (t) . 
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Fig. 3 Trajectory model of atom gravimeters: 
a - fountain type, b - release type. 
Also shown are weighting functions by acceleration 
w g (t) and the disturbance caused by the constant 
gradient Ag~,(t). 



which according to (24 1 and ( 25 1 translates into the 



following measured gravity: 

( 7 

9 = 9a ± 7 I z i + VqT ± Y2 9 ° T ' 



(28) 



This agrees with the result obtained also in 10 15 . 



At the effective measurement height h e s the measured 
gravity equals actual gravity: 



3o ± 7 h e ff- 



(29) 



Comparison of ( 29 1 and ( 28 ) leads to the conclusion 



that the effective measurement height of atom gravime- 
ters is 



7 2 

*cff =zi + V Q T ± ^9qT 



(30) 



For the release gravimeters, let h be the length of the 
idle part of the trajectory, and H be the total trajectory 
length. Obviously, h = z\, H = Z3. As Vo — gozi, 
T=(t 3 - h)/2, t\ = 2 Zl /g , t\ = 2z 3 /g , we get 



h eff = —(H + h) + —VHh. 
24 V ' 12 



(31) 



For the fountain gravimeters, the initial velocity Vq 
equals goT, as the atoms reach the apex in time T. 
The total trajectory height H = g T 2 /2 leads to the 
following location of h e g above the initial position: 



Example 1. Constant vertical gravity gradient and the 
effective measurement height 

We consider two realizations of atom gravimeters. In 
the fountain gravimeters the measurement takes place 
on both upward and downward parts of the trajectory, 
each part has the duration of T (figj3^,). In the release 
gravimeters the measurement takes place only on the 
downward part of the trajectory during the time inter- 
val 2T (fig|3]o) . For both types we put the coordinate 
origin in the atoms' resting position and co-direct z- 
axis with the atoms' initial velocity. The measurement 
in both cases starts at the moment t\ when separation 
from the origin is Z\ and the initial velocity is Vo- The 
constant gradient changes gravity with height like 



g(z) = g a ±j z. 



(26) 



From here on, the upper sign in the ± or =p symbols 
correspond to the release gravimeters, while the lower 
sign corresponds to the fountain gravimeters. In the 



time domain the ( 26 ) is 



g(t) = 3o ± 7 (z 1 + V Q t±g t 2 /2), 



(27) 



Kft = 7.H + z 1 . 
6 



(32) 



As the gravity at the apex of the trajectory is go — 
j(H + zi), the same point is at the following distance 
below the apex: 



h' 



eff 



-H. 



(33) 



The gravity at the effective measurement height corre- 
sponds to the measurement with no vertical gradient 
correction. 



Example 2. Linear vertical gravity gradient 

At some gravimetric sites the vertical gravity gradi- 
ent varies significantly over the free-fall trajectory and 
cannot be considered a constant. The linear gradient 
7i + 72 z changes gravity like 



g( z ) = 9a ± (71 +72z)^- 



(34) 



Unlike the constant gradient, the parameters 71 and 72 
depend on the coordinate system used to analyze the 



correction. Expressing ( 34 1 in terms of time gives us 
g(t) = g + A + Ait + A 2 t 



A,t 3 



AA\ 



(35) 
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Fig. 4 Disturbance caused by the self-attraction 16 and its 
approximation by the 6-th degree polynomial (upper plot), 
and the residuals of the approximation (lower plot). 



where 

A = ±( 7l z! + j 2 zf), Ai = ±(7iVb + 2j 2 V z 1 ) 
M = 7190/2 + 7 2ffo-zi ± IiVq 
A-a = 7250 V Q t 3 , A 4 = ±725o/4. 



Applying the moments (251 to ( 35 1 we obtain the effect 



of the linear gravity gradient as 



Ag = A Q + A\T - 



\a 2 t 2 

6 



A,T 3 



31 
15" 



AaT 



(36) 



Example 3. Self-attraction of atom gravimeter 

The self-attraction of absolute gravimeter can produce 
a complex-shaped disturbance, for which the approxi- 
mations by low-degree polynomials, like in the previ- 
ous examples, are not sufficient. The figure [3] shows the 
self-attraction of the atom gravimeter [16] along the 
atoms' trajectory and its approximation by the 6-th 
degree polynomial. To evaluate the effect of this dis- 
turbance we used moments C n up to the n=12. The 
calculated effect for the inital velocity Vb=0.25 ms _1 
was —1.27 /LtGal, which agrees with the result obtained 
in 16 by another method. The uncertainty introduced 



by the approximation can be estimated with formula 
( 10 ) by treating the residuals e(t) (figQ as disturbance: 



r-2T /-2T 

e(t)w g (t)dt <m&x\e(t)\ w g (t)dt = max\e(t)\. (37) 

D JO 



The result holds true for any type of absolute gravime- 
ter: the uncertainty of the correction caused by the ap- 
proximation of the gravity disturbance does not exceed 
the maximum error of the approximation. In our case 
(figQ the approximation led to the error of less than 
0.05 //Gal. For atom gravimeters the approximation can 
be performed so that the correction uncertainty would 
not exceed max|e(i)|/4. The detailed analysis of this 
issue is out of the scope of the present publication and 
will be presented elsewhere. 



4 Response functions of atom gravimeter 
implementing Bloch oscillations 



In this type of gravimeters 17 the atoms are thrown up 
vertically and travel to the apex of the trajectory during 
the time T, where they hover for the time X" tossed 
by the Bloch oscillations, and then drop back to the 
original position. The formula ([3| for this instrument 
becomes 



17 



02-03 + 04)/[fcT(r + T / )]. 



(38) 



As per section [2j this gravimeter can be modeled by 
three double integrators with integration times of T, 
T + T" , and 2T + T". Impulse response functions of 
this gravimeter are shown on the fig[2jD. As T < T', 
the gravimeter performs near-uniform averaging of the 



17 



atoms' acceleration. The instrument described in 
has the trajectory length of only 0.8 mm, so the correc- 
tions for the constant and the linear gravity gradients 
are minuscule, and the correction for the self-attraction 
can be taken as the disturbance in any point of the 
trajectory taken with the opposite sign. 

The gravimeter's frequency response is Fourier trans- 
form of its impulse response function 18 . Due to its 



almost-uniform impulse response, the atom gravime- 
ter with Bloch oscillations is equivalent to the first- 
order low-pass filter. By comparison (figj5|, the 3-pulse 
atom gravimeter is equivalent to the second order, while 
the corner-cube gravimeter is equivalent to the third- 
order low-pass filter [l8j (fig[5|. These characteristics 
follow from the gravimeters' logic in translating the 
input gravity to the measurement result and do not 
include any additional vibration shielding that instru- 
ments may possess. 

5 Conclusions 

We analyzed atom gravimeters as LTI systems, found 
their impulse response functions and applied them to 
evaluation of certain disturbances. The following con- 
clusions sum up the analysis. 
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Fig. 5 Amplitude of the frequency response function of the 
atom gravimeter with Bloch oscillations (upper plot) com- 
pared to the regular 3-pulse atom gravimeter (middle plot). 
Lower plot: corner-cube gravimeter with multiple levels. Total 
measurement interval: 0.16 s. 



1. Atom gravimeters are equivalent to non-multi- level 
corner-cube gravimeters in translating the gravity 
signal into the measurement result. 

2. Weighting function of atom gravimeter by accel- 
eration can be determined as impulse response of 
the LTI system consisting of several double integra- 
tors. The acceleration, velocity, and displacement 
weighting functions are successive derivatives chang- 
ing sign on every succession. 

3. The effect of a disturbance on atom gravimeter can 
be found by replacing the time powers in the dis- 
turbance expansion with corresponding moments of 
the weighting function. 

4. The effective measurement height of atom gravime- 
ters is located below the apex of the trajectory on 
1/6 of its total length for the fountain, and on about 
7/24 ditto for the release type gravimeters. 

5. Error in the analysis of a disturbance arising from 
the disturbance approximation does not exceed the 
maximum error of the approximation. 

6. With respect to the vibration disturbances the atom 
gravimeter with Bloch oscillations is equivalent to 
the first-order low-pass filter, while the 3-pulse atom 
gravimeter is equivalent to the second order low- 
pass filter. 
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